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Abstract 

We consider the growth of cosmological perturbations to the energy density of dark matter 
during matter domination when dark matter is a scalar field that has undergone Bose-Einstein 
condensation. We study these inhomogeneities within the framework of both Newtonian gravity, 
where the calculation and results are more transparent, and General Relativity. The direction 
we take is to derive analytical expressions, which can be obtained in the small pressure limit. 
Throughout we compare our results to those of the standard cosmology, where dark matter is 
assumed pressureless, using our analytical expressions to showcase precise differences. We find, 
compared to the standard cosmology, that Bose-Einstein condensate dark matter leads to a scale 
factor, gravitational potential and density contrast that increase at faster rates. 



1 Introduction 



Astrophysical observations indicate that 23% of the energy density of the Universe is of an unknown 
nonbaryonic form, known as dark matter. While a precise explanation remains elusive, a weakly 
interacting, nonrelativistic massive particle is favored. Leading candidates include super symmetric 
particles, such as neutralinos, and the axion, originally proposed to solve the strong CP problem 
in QCD. 

Dark matter is often modeled as a pressureless, nonrelativistic particle, known as cold dark 
matter. While this model has achieved significant success, in particular with early universe and 
large-scale cosmology, it meets with difficulty on galactic scales. Cold dark matter simulations of 
galactic halo formation predict density profiles with a central cusp [l|, while observations indicate 
constant density cores [2j . Scalar field dark matter that has undergone Bose-Einstein condensation 
[3 -10 ] has been considered as a solution to this problem since the resulting density profiles agree 



with observed rotation curves (5L fllTllli] . 

At low temperatures, where wave aspects dominate, a many-body system of bosons exhibits 
Bose enhancement, whereby bosons favor joining highly populated states. As a result, bosons 
can pile up in the same ground state forming a coherent matter wave of macroscopic size known 
as a Bose-Einstein condensate (BEC). This will occur when the thermal de Broglie wavelength, 
A = v /2irti 2 /2mk B T, be gins to exceed the interparticle spacing, n 1//3 , so that wave functions of in- 
dividual bosons begin to overlap, where m is the boson mass and n its number density. Equivalently, 
Bose-Einstein condensation occurs when the temperature, T, drops below the critical temperature, 

T - I ^ n 2/3 fl) 

where C( x ) * s the Riemann-Zeta function. 

The experimental realization of trapped BECs in dilute alkaline atoms in 1995 has led to 
a renewed interest in BECs, which has been a subject that unifies many disciplines, for example 
neutron stars 0, superconductivity Q and, what is our interest here, dark matter |3l4l0[|. A 
great virtue in the study of BECs is that in ultracold atomic systems important system parameters, 
including the interaction strength between particles and the dimensionality of the system, can be 
tuned precisely. This opens up the possibility of using ultracold atomic systems in the laboratory 
to simulate phenomena on galactic and cosmological scales. Possible examples include creating a 
system where the BEC is subject to so-called electromagnetically induced "gravity" with a 1/r 



interatomic attractive potential [17|] and creating a controlled explosion of atoms by suddenly 
making the s-wave scattering length negative, a phenomenon dubbed "bosenova" because of its 
resemblance, on a vastly lower energy scale, to the core collapse in a supernova fl8| . 

Our focus here is with the cosmological applications of BECs, in particular BEC dark matter. 
In addition to the density profile and rotation curves mentioned above, investigations of BEC dark 
matter include the study of vortex formation [f], 19-2 2|. a dditional aspects of galactic structure (23| . 
Bose-Einstein condensation in the early universe [24-26] and axions |27]. Recently a study of its 
cosmology was initiated by Harko [28| and Chavanis (29I] (see also 25|, 30-32]). They investigated the 
evolution of inhomogeneities in the dark matter energy density. Such inhomogeneities eventually 
become nonlinear and lead to galaxy formation and indirectly affect anisotropies in the cosmic 
microwave background radiation. They derived evolution equations for the density contrast in 
Newtonian gravity [2^] and post-Newtonian gravity 0, [3] and presented numerical solutions to 
these equations. 

In this work we also analyze the evolution of cosmological inhomogeneities after dark matter 
has undergone Bose-Einstein condensation. Our departure from previous results is twofold. First, 
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we present simple analytical solutions which allow for a precise understanding of how BEC dark 
matter differs from standard cold dark matter. Since dark matter is believed to have a small 
pressure (indeed, it is often modeled as having zero pressure), by taking the small pressure limit 
such analytical solutions are obtainable. Second, while we will begin with Newtonian gravity, we 
derive the evolution equations using the complete theory of General Relativity. The use of General 
Relativity is necessary when considering superhorizon perturbations, which are beyond the reach 
of Newtonian gravity It is also necessary for considering anisotropies in the radiation spectrum. 
We do not study anisotropies here, as it is outside the scope of our work, but such a study would 
be interesting and important and would rely on our results. 

In the next section we review the Gross-Pitaevskii equation coupled to the Poisson equation 
and the Thomas-Fermi approximation, which has become the standard framework for describing 
a gravitating BEC. We then derive the equation of state for BEC dark matter and use it to 
analyze the homogeneous, unperturbed cosmology. In section [3] we derive the evolution equations 
for inhomogeneities during matter domination, first in Newtonian gravity and then in General 
Relativity. Throughout we compare our results to those of the standard cosmology, where dark 
matter is pressureless. When making these comparisons we shall refer to such dark matter as 
standard cold dark matter (SCDM). We conclude in section [U 



2 Bose-Einstein Condensate Dark Matter 
2.1 Hydrodynamic Description 

We assume that dark matter is composed of scalar bosons, of mass m, that have undergone a 
phase transition toward Bose-Einstein condensation. To describe the BEC we employ the standard 
symmetry-breaking mean field approach, which is expected to be valid for systems with a sufficiently 
large number of particles and at temperatures far below the BEC transition temperature. In this 
approach we may start with the Gross-Pitaevskii energy functional: 



E\M= / drx 



^-\V^(t,x)\ 2 + ^Vo\iP(t,x)\ 4 + ~mV G (t,x)m,x)f 



(2) 



where ip(t, x) is the order parameter, or macroscopic wave function, describing the BEC and is 
normalized such that \ip\ 2 is the number density. The first term is the standard kinetic energy term 
of nonrelativistic quantum mechanics. The second term represents a quartic, contact interaction 
with strength 

A-nh 2 a s , . 

V = -, 3 

m 

where a s is the s-wave scattering length, which we take to be positive (a s > 0). The third term is 
the gravitational potential, 

V G (t,x) = -Gm f dV 1 ^'^' 2 , (4) 
J \x — x'\ z 

which satisfies Poisson's equation: 

V 2 V G (t,x) = 4-KGm\tp(t, x)\ 2 . (5) 

By including the gravitational potential, ([2]) describes a BEC coupled to gravity. 

In the mean field approach one ignores high order correlations due to bosonic quantum field 
fluctuations. This allows the BEC to be described by the equations of motion that follow from 
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variation of ([2]), under the constraint that the total number of particles is conserved: 



Bib ft 2 

ih^ = -—V 2 ^ + V \^ + rnV G ^-^, (6) 

where p is the chemical potential. This equation is known as the Gross-Pitaevskii equation. It 
may be written in the hydrodynamic representation, which is more useful for our purposes, by 
separating the wave function into its modulus and phase, 

4,(t,x) = m,x)\e iS ^\ (7) 

both of which are real, and then describing the condensate in terms of its energy density and local 
velocity: 

p(t,x) = mc 2 \ip(t,x)\ 2 , v(t,x) = — VS(t,x). (8) 

m 

In terms of these variables the Gross-Pitaevskii equation © and the Poisson equation ([5]) become 

% = -V • ( P v) (9a) 

-at = {-p v n + 2 V{V } + ^ 

W G = ^p. (9c) 

Aside from the first term on the right hand side of (|9b|) . the top two equations comprise the 
hydrodynamic description of the condensate since (|9a|) is the continuity equation and (|9b|) is the 
Euler equation from classical fluid dynamics. The first term on the right hand side of (|9b|) is 
called the quantum pressure term. Unfortunately this term often makes analytic solutions difficult 
to come by. It may be traced to the first term in ([2]), which originates from the uncertainty 
principal and hence cannot find its analog in classical physics. Because this term contains the 
gradient of the energy density, as the number of particles in the condensate, or equivalently the 
size of the wave function, increases, the quantum kinetic energy becomes negligible compared to 
other energy contributions except near boundaries of the condensate. Neglecting the quantum 
pressure term is known as the Thomas-Fermi approximation (for a more precise definition of the 



V - V P + oV(^) + ^Vp + VU G (9b) 



gravitational Thomas-Fermi regime see [2l|])- It is a common practice to employ the Thomas- Fermi 
approximation in the study of perturbations to BEC densities something we shall do in the 
next section. 



2.2 Equation of State 

The equation of state for a fluid relates the pressure to the energy density, p = p(p), and is of 
fundamental importance in cosmology. It may be obtained for a homogeneous BEC by ignoring 
the gravitational potential in ([2]) and taking ip(t, r) = ipo to be real and constant, where ib 2 = N/V 
is the number density. From ([2]) the BEC has energy 



E 



I N 2 



(10) 



The pressure is then 



dE 
dV 



N 



(11) 
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In a BEC each particle contributes an energy mc 2 , so the energy density is mc 2 times the number 
density, p = mc 2 ipQ, and we find the equation of state 

V 2 2irh 2 a s 2 _ 2 
P= 2ri& P =^W P = Xp > (12) 

where A = 2irh 2 a s /m s c . 

Dark matter is thought to be cold and nearly pressureless. The standard assumption is that it 
is a pressureless, perfect fluid, which we will refer to as standard cold dark matter (SCDM), with 
equation of state Pscdm = 0. BEC dark matter has nonzero pressure and the nontrivial equation 
of state (|12p. We will analyze BEC dark matter using the dimensionless quantity 

w = I = \ p . (13) 
P 

When solving for approximate, analytical solutions, we will take w to be a small perturbation 
around the SCDM solution u^scdm = 0. The SCDM results can be obtained by setting w = 0. 

In general w is not constant, but we will make use of it evaluated at its (constant) present-day 
value ivq. Further, we can introduce the dark matter fraction f^DM = P/pci where p c is the critical 
energy density necessary for a flat universe [jjji]. Then 

w = A^dm,oPc,o (14) 
where Odm,o and p Ct o are the present-day values. 

2.3 Homogeneous, Isotropic Cosmology 

If we ignore perturbations, the Universe is well known to be flat, isotropic and homogeneous on the 
distance scales of interest [34]. It may be described by the evolution of the scale factor, a(t), which 



evolves according to the Friedmann equations 33] 



where H = a /a is the Hubble parameter and a dot denotes a time derivative. Here, and from 
now on, the subscripted on the energy density, po, and the pressure, po, indicate that these 
are unperturbed, background quantities. The cosmology of BEC dark matter follows from the 
Friedmann equations and the equation of state. In section 13.11 when studying inhomogeneities in 
Newtonian gravity, we will derive the Friedmann equations directly from (J9j) . Here we simply quote 
their well known form. 

Our goal in this subsection is to determine the evolution of the scale factor during matter 
domination when BEC dark matter, with equation of state f)13f) . dominates the total energy density 
of the Universe and perturbations have been ignored (perturbations will be considered in the next 
section). Using (fT3l) and the second equation in (fT5l) we have (28| 

A 

where A, since it is the exponential of an arbitrary constant, is positive, but otherwise arbitrary. 
It may be fixed by requiring the energy density to have its present-day value, po = po,o ; when the 
scale factor has its present-day value, a = clq, leading to A = po^a^/il + Apo.o)- With this we may 

/ s Po,o(l ~ W ) 
P0ia) = (a/aor-Wo> (17) 



Po( a ) = -3 — rr> ( 16 ) 



rewrite (1161) as 
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where for convenience we denned Wo = wq/ [1 + wo) . Using the solution (fT7|) and the first equation 
in (HH) we obtain Q 



n DM ,o(l - W )H (t - t') 



(y - tan 1 y) , 



(18) 



where t' is an arbitrary constant, and 




(19) 



t' may be fixed by applying an initial condition. Initially, in the very early Universe when the 
temperature was sufficiently high, BEC dark matter will not yet have condensed. As the Universe 
expands and cools, dark matter will eventually begin to condense, taking a finite period of time 
to complete [26|. Once complete, (fT7|) is valid, so that completion occurs during a > aoW^ 3 . 
The initial condition we apply, which was used in 28, is a(t = 0) = aoW^ 3 , or equivalently 
y (t = 0) = 0, which sets t' = 0. 

So far everything we have done is exact. We will now use the approximation wq <C 1, which is 
the statement that the present-day pressure is small. It follows from this that W = wo + 0{wq) 
and y ~ yj {a/ao) 3 /wo is large. Expanding around large y and small wo we obtain 



v /^DM,o(^oi) = -yw + O(wl) 



71" 1 , 

2 y 



This equation may be solved for a and then expanded around small wo to obtain 

2vr 24 + vr 2 



ait) _ ^90dm,o x1/3 



ao 



(H t) 2/3 + 



(20) 



+ 0(™f). (21) 



This is the desired equation expressing the evolution of the scale factor in the small pressure limit 
during matter domination for BEC dark matter. We note that upon setting wo = we immediately 
obtain the SCDM result: oscdm oc t 2 / 3 . The additional terms are modifications taking into account 
the nonzero pressure of BEC dark matter. 

Using the WMAP result f2DM,o = 0.228 we have plotted the scale factor in figured! The 
dashed line on the bottom is the SCDM solution. The solid lines are our approximate, analytical 
solutions (f2Tj) for various wo (see caption) and the dotted lines the exact solutions (fT8|) for the 
same wo. As can be seen, as Hot gets smaller the analytical solutions become less accurate. The 
reason for this can be traced to our assumption that y in (|19p becomes large for small wo- As a/ao 
decreases, this assumption becomes less valid. In (12ip this manifests itself as the magnitude of 
the wo terms increasing for smaller Hot, making the truncated expansion in (|21|) less valid. Since 
present-day has been defined as a/ao = 1, we can see that BEC dark matter leads to a larger scale 
factor than for SCDM. 



3 Inhomogeneities 

In section [231 we considered the homogeneous, isotropic Universe. In this section we perturb around 
this Universe with the goal of determining the evolution of the perturbations, or inhomogeneities, to 
the energy density of BEC dark matter. We begin in the following subsection with inhomogeneities 
in Newtonian gravity For nonrelativistic matter and subhorizon perturbations, Newtonian gravity 
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Figure 1: The evolution of the scale factor is plotted during an epoch of matter domination. The 
dashed line corresponds to the SCDM result. The solid lines are our approximate, analytical results 
(|2ip for BEC dark matter. From top to bottom they correspond to wo = 0.1, 0.01 and 0.001. The 
dotted lines are the exact results (fTHj) for the same values of wq. 



is (nearly) sufficient for determining the leading order solutions [34] . In section 13.21 we solve for 
the evolution of inhomogeneities in General Relativity. While the use of the complete gravitational 
theory of General Relativity has the benefit of verifying our Newtonian solutions, it also intro- 
duces relativistic corrections and allows for the consideration of superhorizon perturbations, which 
Newtonian gravity is incapable of describing. While we could present only the General Relativistic 
results, we find the Newtonian analysis more transparent and some discussions in section 13.11 are 
necessary for justifying a straight forward application of General Relativity in section [3~2l 

3.1 Inhomogeneities in Newtonian Gravity 

In the absence of perturbations, the Universe is well known to be homogeneous and isotropic on 



the distance scales of interest 34j]. This means the energy density of matter does not vary over 



space and the velocity of matter is due only to expansion and obeys the Hubble law: 

Po = Po(t), v = v (t, x) = H(t)x, (22) 

where the subscripted 0's indicate that these are background, unperturbed quantities. Before 
perturbing them, there are two important equations we can derive. Plugging these quantities into 
the continuity equation ([9a]) and the divergence of the Euler equation (f9b|) combined with the 
Poisson equation (|9c|) we find 

Po = -3tfp , H + H 2 = ~^-po, (23) 

which are the Friedmann equations and, as promised in section [2T2j we have derived them directly 
from 0. Note that compared to (|15p . the Friedmann equations here are missing the pressure 
term. This is because these are the nonrelativistic Friedmann equations, derived within Newtonian 
gravity. If we were to write our equations in terms of the mass density instead of the energy density, 
then, as can be seen from (|15p . only the pressure would contain factors of c, showing that it is a 
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relativistic correction. An alternative form for these equations will soon be of use. Below we will 
introduce the scale factor, a(t), which will be related to the Hubble parameter by H = a/a, where a 
dot denotes a time derivative. In terms of the scale factor the two equations in (j23[) can be written 



2 8irGa 2 AirGa , d . 

a =—^~Po, a = — ^Po- (24) 

We now perturb the above quantities around their background values: 

p(t,x) = p (t) + 5p(t,x) (25a) 

v(t, x) = v {t, x) + Sv(t, x) = H(t)x + 5v(t, x) (25b) 

V G (t, x) = V G0 (t, x) + 5V G0 (t, x). (25c) 

In terms of these perturbed quantities, Q becomes 

(5p)x = ~P(Nx -dv-Vx- (Spv ) (26a) 

Si = -^-l^ x (yls p ) + v x (v • 8v) + — VJp + V X 5V G (26b) 
2m po po 

Vl5V G = ^f5p, (26c) 

where we've used the unperturbed version of ([9]) to cancel terms and w was defined in (|13p . Soon 
we will Fourier transform these equations to facilitate solving them. However, in their present form 
they will mix Fourier modes. To avoid this we move to comoving coordinates, q, given by x = a(t)q, 
which requires transforming derivatives as 34] 

d \ ( d \ 1 

al) -^o-V*, V x = -V q . (27) 



Then (|261) becomes 



dt I \dt 



S = —V-SH (28a) 



a 

-5v = - — ^V(V 2 5) + H5v + ^^V5 + - V5V G (28b) 
2m a 6 a a 

V 2 5V G = ^f-a 2 p 5, (28c) 



where 



8 = *± (29) 
Po 



is the density contrast. In (|28p we refrained from writing subscripted q's and, for the Euler equation 
(|28bp . used the fact that the velocity of a BEC is irrotational, as follows from ([8]). These equations 
may be combined by taking the divergence of the Euler equation (|28b|) and then subbing into it 
the continuity equation (I28ap and the Poisson equation (I28c|) . giving 

h 2 1 . ■ 2wc 2 9 4vrG , . 

V 4 5 - 2H8 + —^V 2 8 + —pot- (30) 



2m o 4 a 



This equation describes the evolution of the density contrast, S, in an expanding universe. Our inter- 
est is to solve this equation and look for growing modes representing the growth of inhomogeneities. 
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Such inhomogeneities will eventually become nonlinear (5 > 1) leading to galaxy formation. We 
focus here on the linear regime (6 < 1), where a perturbative analysis is accurate, paying particular 
attention to the rate at which inhomogeneities in BEC dark matter grow compared to SCDM. 
To solve (p0|) we begin by Fourier transforming the density contrast, 



5(t, x) 



so that (j30j) becomes 



5 k + 2H5 k + 



2m a 4 



d A k 
(2tt) 2 / 3 



+ 2wc 2 



5 k {t)e 



ik-x 



4ttG 

— —po ) Ok 



0. 



(31) 



(32) 



At present, to solve (|32j) one needs to specify the exact time dependence of the scale factor a(t). 
To avoid this, we can transform the independent variable from cosmic time, t, to the scale factor, 
a, with the help of pi]) , to obtain [29]] 



d 2 5 k 3 d5 k 



+ 



(. 



2 1.4 



h z k 



+ 



wc k 



.41.2 



da 2 2a da 2a? \8irGm 2 a 4 po 2nGa 2 pi 



1 5, 



0. 



(33) 



where now 5 k (a), po(«) and if (a) are functions of the scale factor. 

For large enough k, the solutions are oscillating sound waves. From (I32p we can see that this 
corresponds to the gravity term being negligible and pressure dominating. For small enough k 
the gravity term dominates and the solutions are growing or decaying. Since we are interested 
in growing solutions, our focus is then on smaller modes. Ideally we would solve (|33p exactly. 
Unfortunately, it does not appear possible to find exact solutions [29]. The problem is due to 
the A; 4 term, which comes from the quantum pressure term in ()9b[) . At the end of section [2TT1 we 
mentioned that the quantum pressure term often causes difficulty for finding analytical solutions 
and that neglecting it is known as the Thomas-Fermi approximation. We will now employ this 
approximation. Note that making the Thomas-Fermi approximation leads to the standard classical 
hydrodynamic equations 34]. We will make heavy use of this fact in the next subsection. 

Making the Thomas-Fermi approximation, (|33p can be written 



d 2 5 k 3 ddk 3 
da 2 2a da 2a? 



ao\ 2 ( k_ 
aJ\ kj 



0, 



(34) 



where ao is the present-day value of the scale factor and kj = aftirGpo/wc 4 ) 1 / 2 is the Jeans mode 
with present day value 



kj, 



o — a o 



2<irGp 



a 



1/2 



(35) 



\ 4c 2 wo 

The Jeans mode gives the exact point of crossover from oscillating solutions to growing and decaying 
solutions. Oscillating solutions occur for modes k > kj while growing and decaying solutions occur 
for modes k < kj. Since kj^ is independent of a, all de pen dence on the scale factor (other than in 



5k) has been explicitly written in (|34p . The solution is 29j] 



Sk(a) = C' kl (^) 1/i J_ 



5/4 



oq\ k_ 
a ) kjn 



a 



a \i/4 



J- 



5/4 



V 2 U ) kjn 



(36) 



where C' kl and C' k2 are arbitrary constants and J±5/4 are Bessel functions of the first kind. In 
figure |2f a) we have dropped the decaying solution (C' k2 = 0) and plotted the growing solution for 
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0.0 0.2 0.4 0.6 0.8 1.0 

a/ao 

Figure 2: In (a) the exact solution (j36h is shown for k = kj$. In (b) the dashed line is the 
wq = SCDM solution, the solid lines are our approximate, analytical solutions ([37|) for, from 
top to bottom, k 2 /kj Q = 1/10, 1/20 and 1/50. The dotted line is the exact solution (|36p for 
1/20. 



k = kjfi. For a < ao, the solution is oscillating, since pressure is dominating over gravity, as can 
be seen from (I34D . For a > ao, gravity begins to dominate and the growing solution emerges. 

As mentioned previously, the direction we take in this paper is to find analytical solutions in 
the small pressure limit that, upon setting wq = 0, reproduce SCDM results. Such a solution can 
be obtained by expanding ([36]) . but we must be careful to expand around a small quantity. Since 
growing solutions occur for k < kj, we will focus on modes that satisfy k <C fcj,o and k/kj t o will be 
our small quantity. For these modes, we find 



5 k (a) = C, 



fcl 



a \ 3k 2 ( a 
— + w ^- — 
a J 2k 2 \ a 



+ Cfc2 



-3/2 



W 



6k \ao 



-7/2 



+ 0(w 2 ), (37) 



where k = wok 2 (and is independent of wq), Cki is an arbitrary constant proportional to C' kl and 
likewise for and C' k2 - One can see clearly that for wq = this equation reproduces the SCDM 
solution: 4, SC DM = C' k \a + Cj^cT 3 / 2 0. 

Dropping the decaying solution (C k 2 = 0) we have plotted the growing solution in figure E^b). 
The dashed line is the wq = SCDM solution, the dotted line is the exact solution (|36|) and the 
solid lines are our approximate, analytical solutions (f37l) for various values of k 2 /k 2 (see caption). 
Here again we find that the approximate solutions become less valid as a/ao decreases. The reason 
is the same as before: as a/ao decreases, the wo term in ([37]) increases, lessening the validity of the 
truncated expansion. We also find that inhomogeneities are larger with BEC dark matter. In a 
universe with BEC dark matter, then, galaxy formation is expected to happen sooner than in the 
ACDM universe 



To write the solution (|37p in terms of cosmic time, t, requires knowledge of how the scale factor 
evolves with time. Since the evolution of the scale factor changes as different components dominate 
the total energy density, we must further specify a particular epoch of the Universe during which 
we wish to determine a(t). In (|2ip we found a{t) during matter domination when BEC dark matter 
dominates the total energy density. It is customary to drop the purely decaying C&2 solution so 
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that, upon subbing ([2T|) into (f3T|) . we have 



s k (t) 



9a/^dm,o 




(#0*)" 



-1/3 



ot"DM,0 



(38) 



Again, we find the SCDM solution when setting wq = 0. 

In (|38p the leading modification to the SCDM solution is positive. Thus, we see analytically that 
BEC dark matter leads to an increased rate for the growth of inhomogeneities compared to SCDM. 
This result is in line with 28, H^], however in those papers the authors determined their complete 
solutions numerically. Here we have obtained analytical solutions for growing modes during matter 
domination. It is well known that the growth of inhomogeneities occurs at an appreciable rate only 
during matter domination 34J], and thus we have focused on this epoch. 

In this subsection we have made a nonrelativistic, Newtonian analysis for the evolution of 
inhomogeneities. However, in our final equation (j38f) we used the relativistic result (|2ip which 
followed from the fully relativistic Friedmann equations (|15p and not the Newtonian Friedmann 
equations (|23j) . Using the Newtonian Friedmann equations would remove the y/wo term. The 
importance of this term, then, begs the question of how important relativistic corrections are for 
the density fraction. One possibility would be to make a post-Newtonian analysis, as was done in 
2~8l . 29]. Post-Newtonian gravity 37, 38] includes the leading relativistic corrections from General 
Relativity, and not just those in the Friedmann equations. We opt instead to make a fully relativistic 
analysis using General Relativity. 



3.2 Inhomogeneities in General Relativity 

From this point forward we set c = 1. We consider only a flat universe and scalar perturbations 
and write the metric in conformal Newtonian gauge as 



ds 2 = a 2 (r]) [(1 + 2^)drj 2 - (1 - 2<f>)5 ij dx i dx j ] 



(39) 



where and $ are scalar perturbations, rj is conformal time related to cosmic time via dt = adrj and 
the x l are comoving coordinates. The Hubble parameter, in terms of conformal time, is given by 
H = a' /a, where a prime will denote differentiation with respect to rj, and the Friedmann equations 
are 



H 2 = 

which may also be written as 

l\ 2 

a 



8ttG 



a 2 po, 



AttG 



a 2 (po + 3p ) , 



IttG n a" 

-5— a Po, — 
6 a 



4ttG 



a 2 (p - 3po). 



(40) 



(41) 



The Hubble length, or horizon, is given by the physical distance H~ l = (aH)^ 1 ~ arj. A 
perturbation with comoving mode k ~ 1/A has a physical wavelength of roughly aX. Thus hrj < 1 
corresponds to superhorizon modes with physical wavelengths longer than the Hubble horizon and 
krj > 1 corresponds to subhorizon modes with physical wavelengths shorter than the Hubble horizon. 
In the previous subsection we considered only subhorizon modes, since this is all a Newtonian 
analysis can accommodate. In this subsection we consider both subhorizon and superhorizon modes. 
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The metric (I39p obeys the Einstein field equations 



Gl = SvtGT;, (42) 

where the Einstein tensor G 1 ^ is a function of the metric and T£ is the stress-energy tensor. The 
perturbed Einstein field equations, 5G^ = 8ttG5T^, are 34J 



4TrGa 2 5T° = V 2 $ - 3H + %<$>) (43a) 
4TrGa 2 5T° i = (y' + Wit) (43b) 



4irGa 2 5T) = i(* - tf)^ 



tf" + ft(2tf' + + (2W 7 + % 2 )<I> + ^V 2 ($ - *) 



Sij, (43c) 



where we have used the standard notation that a Latin index refers to spatial components only and 
a comma denotes a partial derivative, e.g. \& j = d^/dx l . In this subsection, as above, we restrict 
our attention to the epoch of matter domination, during which BEC dark matter is dominating 
the total energy density of the Universe. Then is the stress-energy tensor for dark matter. 
As explained in the previous subsection, we are able to drop the quantum pressure term in the 
Euler equation (|9b|) by making the Thomas-Fermi approximation. This approximation reduces the 
equations ([9]) to the standard classical hydrodynamic equations of a perfect fluid. Such equations 
may be derived from the stress-energy tensor for a perfect fluid, 

r; = ( P + P ) u \-^, (44) 

where u M is the 4-velocity of the fluid. Conservation of the stress-energy tensor and the nonrel- 
ativistic limit reproduces exactly the continuity equation (J9aj and the Euler equation (|9bj) if one 
drops the quantum pressure term and makes the identification p = Vop 2 /2m 2 c 4 , which is identical 
to the equation of state (112j) . which was derived through other means. Perturbations to the perfect 



fluid stress-energy tensor are given by 34] 

5T° = 5p, 5T* = a(p +p )5v\ 5Tl = -5p5{. (45) 

Since the perturbations to the stress-energy tensor are diagonal, (|43cp tells us <I> = ^> . In the 
following, then, we shall label all scalar perturbations with <I>. 

BEC dark matter satisfies the equation of state (fl~3j) . so that 5p = 2w5p. Combining (|43ap and 
(|43bj) we obtain 

$1 + 3(1 + 2w)n<5>' k + \TH! + n 2 + 2w{m 2 + k 2 )]<5> k = 0, (46) 
and from (|43ap alone 



(47) 



where we made use of the Friedmann equations, 6% = 5pk/po is the density contrast and we have 
Fourier transformed both <1> and 5. If we are able to solve (|46p for <£fc, then we may place it into 
(I47p to obtain 5k- Note that, in general, the 5k in (I47p is the density contrast for the perturbation 
to the total energy density. Since we are restricting ourselves to the epoch of matter domination, 
the total energy density is the energy density of dark matter and the 5k in (|47|) is the same 5k that 
we solved for in the previous subsection, the inhomogeneity to BEC dark matter during matter 
domination. 
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To solve (|46p we first transform the independent variable from conformal time, 77, to the scale 
factor, a, yielding 

d 2 $ k 1 . _ d$ k . 3w r i + k 2 N| $ 



da 2 



+ 



47rGa 2 p / 



0. 



Using the results in section 

52 



we can exchange w for its present-day value, wq, and obtain 



H 

da 2 4a 



14 + wo 



30 



9a 



1 a fc 2 \ 
1 + 2^p)^ + °^ ) 



0. 



(48) 



(49) 



a-a («/ a o) 5 

We could have written this equation in its entirety, and not just through order wq. However, it 
does not appear possible to solve the complete equation analytically. Instead we solve the equation 
in the small pressure limit. This allows use to solve for the solution expanded around small wq. 
We find 




-9/2 



21 

22 



-11/2' 



(50) 



+ OK 2 ) 



where CL and C' k2 are arbitrary constants. This equation gives the evolution of the gravitational 
potential for any a. During matter domination we may write it in terms of cosmic time by using 
(|2"Tjh It is customary to drop the purely decaying C k 2 solution so that we have 



$ & = C' k J 1 + w 



(— 



2/3 



3fe 2 , N 
2k 2 V ' 



■4/3 



1/3 



\9^DM,0/ 



(Hat 



-2 



+ 0(^ /2 ). (51) 



During radiation domination, which preceded matter domination, <]?£, is constant for super- 
horizon modes 39]. Since the Hubble horizon expands faster than the physical wavelength of a 
perturbation, the wavelength of the perturbation eventually becomes subhorizon. If this occurs 
during radiation domination, the perturbation decays quickly to zero [39]. If instead the perturba- 
tion survives into matter domination, then it evolves according to (15 1 f) . We find immediately, upon 
setting wq = 0, that (|5ip reproduces the SCDM solution of a constant gravitational potential for 
both superhorizon and subhorizon modes. BEC dark matter introduces ^-dependent corrections, 
so that the gravitational potential is no longer constant. 

With the solution (I50p for the gravitational potential we can obtain the density contrast using 
(|47p . To do so we first transform the independent variable in (|47p from conformal time to the scale 
factor and then, analogously to how we arrived at (I49p . we use the results in section [2T3l to find 

-3 N 



-2 ^ ad a <f> k + 



1 + 



k z 



4k 2 ao 




WQ + Wq — 



(52) 



Into this equation we plug (|50j) to obtain the desired result: 

"3^2 / ... x -1 / ... x 2 

2k 2 

-3/2 



ki 



« 



+ Wq 



7 



4/c 2 
k 2 




6k 2 



-7/2 



l + ^o + ^o 



63 
22 



283 




"66" 


( a o. 


") 


-3 


1 









1 + w + 4w 
9/2" 



-3 



(53) 



+ 0(w 2 
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where Cf-i = —k 2 C' kl /2k 2 (l — wq) and likewise for Ck2 and C' k2 . This result may be compared to 
the Newtonian result (J37J) . The fully relativistic solution ([53]) reproduces the Newtonian terms (|37p 
exactly. We may rewrite (|5ip in terms of cosmic time by using (|2ip . It is customary to drop the 
purely decaying Cik solution, so that we have 



$k = Cfci 



/ 90dm,o x 1 



{Hot) 2 " + 



2tt 



wr 1/3 



+ W Q 



k 2 



1/3 



k 2 \ 3n 2 



Hoty 2/3 - I^Vo*r 4/3 



+ 



4p / 



DM,0 / 
4 \ 1/3 



81^DM,0 



(54) 



k 2 \9S7dm 



1 + Wq + W o 
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9^DM,i 



,-2 



During radiation domination, 8f. for SCDM is constant for superhorizon modes and grows at 
most logarithmically for subhorizon modes [3^]. The growth of inhomogeneities is largest during 
matter domination, during which 8^ for BEC dark matter evolves according to (|54p in the small 
pressure limit. For wq = 0, (|54p reproduces the SCDM solution. The dominant additional terms 
are positive and thus the growth rate of inhomogeneities increases for BEC dark matter compared 
to SCDM. For subhorizon modes, kr] ~ kt 1 " 3> 1 and we can drop the k~ 2 terms. This case was 
already discussed in the previous subsection. For superhorizon modes, krj ~ kt l " <C 1 and we can 
drop all but the k~ 2 terms: 



S k = -2C' kl 



1 + 2wo + u; 
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90 D 



M,0 



<H Q t) 



-2 



(superhorizon). 



(55) 



From this solution we can see that during matter domination superhorizon inhomogeneities for 
SCDM are constant, while for BEC dark matter there is growth. 



4 Conclusion 

In this work we studied the resulting cosmology when dark matter is a scalar field that has undergone 
Bose-Einstein condensation. Such a model of dark matter has been shown to be in better agreement 
with the density profiles of galactic halos than standard cold dark matter in the ACDM model. We 
focused on the growth of inhomogeneities, i.e. perturbations to the dark matter energy density. 
Since, as is well known, such perturbations only grow appreciably during matter domination, we 
considered only this epoch. 

Such an analysis had been considered in earlier work 0, [2^| within the context of Newtonian 
and post-Newtonian gravity. Our analysis differed from these papers principally in two ways. 
First, the direction we took was to derive simple, analytical formulas that clearly showcased the 
modifications BEC dark matter produces in the standard cold dark matter solutions. Second, 
while we did make our analysis first within Newtonian gravity, we then used the complete theory 
of General Relativity, finding solutions for the scale factor, gravitational potential and density 
contrast. We found analytically that each one of these quantities increases at a faster rate compared 
to when dark matter is in the form of standard cold matter, consistent with the numerical results 



m 



Our fully relativistic solutions for the gravitational potential in f)51 1) and the density contrast 
in (I53p are valid for both subhorizon and superhorizon perturbations, the latter of which is beyond 
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the reach of Newtonian gravity. These fully relativistic solutions are also necessary for studying 
anisotropies in the radiation spectrum. While we did not include an analysis of anisotropies here, 
such an analysis we expect to lead to new and interesting physics and would rely on our results. 
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